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We outline the recent results on the ground state for a class of one- and two-dimensional frustrated 
quantum spin models with competing ferro(F)- and antiferromagnetic (AF) interactions. Frustrated 
, spin systems are known to have many interesting properties due to large quantum fluctuations. As 

a result of these fluctuations usual mean-field approach gives quite crude (if not false) description 
of these systems. Therefore, exactly solvable models are very instructive at investigations of such 
systems. The exact ground state wave function of proposed models has a structure of the valence- 
bond state (VBS) type. One of the ID model describes the transition line between the F and AF 
phase. The exact singlet ground state on this line has a double-spiral ordering. Using different 
approximation methods we study the magnetization curve in the AF phase. The second considered 
set of the ID and 2D models has an exact non-degenerate ground state with exponentially decaying 
spin correlations. We also proposed the ID and 2D electronic models with exact ground state 
represented in terms of singlet bond functions which are the generalization of the RVB functions 
^ including ionic states. 
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There is currently much interest in quantum spin systems that exhibit frustrationS. This has been stimulated in 
particular by the study of the magnetic properties of the cuprates which become high-Tc superconductors when doped. 
Frustrated spin systems are known to have many interesting properties that are quite unlike those of conventional 
^ magnetic systems. The simplest model of such kind is the Heisenberg_spin chain with nearest- and next-nearest 
O [ neighbor interactions Ji and J2. This model is well studied for Ji, J2 > OcTu. In uarticular, it has been found that at 
J2 = O.24J1 the transition from the gapless^state to the dimerized one takes placcfl. The point J2 — Ji/2 corresponds 
to the well-known Majumdar-Ghosh modelQ for which the exact ground state consists of dimerized singlets and there 
is a gap in the spectrum of excited states. 
^ ' Less studied is frustrated spin models with competing interactions of ferro- and antiferromagnetic types. The 
\^ , physical interest for these models is connected with the study of the real compounds containing CuO chains with 
edge-sharing CUO4 units, like LaQCasCu2404i, Li2Cu02 and (70212(^^50100. In these compounds Cu — O — Cu bond 
angle is nearly 90° and the nearest-neighbor Cu — Cu spin interaction is ferromagnetic according to Goodenough- 
Kanamori- Anderson rulecl while the next-nearest-neighbor interaction is antiferromagnetic. The magnetic properties 
—\ of these models are very different from those for models with pure antiferromagnetic interactions. Their ground 
^ . state can be either ferromagnetic or singlet depending on the relation between the exchange integrals. One of the 
^ ■ most interesting problems related to these models is the character of the transition between ferromagnetic (F) and 

antiferromagnetic (AF) phases. 
^ ' Of a special importance are models for which it is possible to construct an exact ground state. Last years a 
' considerable progress has been achieved by using so-called Matrix-Product (MP) form of the ground state wave 
functioroEj. The ground state wave function in the MP method is represented by Trace of a product of matrices 
Q ' describing single-site states. The MP ground state has a structure of the type where each neighboring pair of spins 
Q has valence bond and, in fact, the MP form is a convenient representation of the valence bond states. Its origin can 
^ ■ be traced to the S=l quantum spin chain with biinear and biquadratic interactionalj. At present, various ID spin 
. ^ models with the exact MP ground state is foundt2l"tHl. 

In this paper we present a set of ID and 2D spin-1/2 models with competing F and AF interactions for which 
H ] the singlet ground state wave function can be found exactly. This function has a special form expressed in terms of 
■ - - 1 auxiliary Bose operators. This form of the wave function is similar to the MP one but with infinite matrices. For 
special values of model parameters it can be reduced to the standard MP form. 

One of models is ID quantum spin model describing the F-AF transition point. Spin correlations in the singlet 
ground state show giant spiral magnetic structure with period equals to the system size. On thej-antiferromagnetic 
side of this point the ground state can be either gapless with incommensurate spiral ground statala or gapped with 
exponential decay of correlatoraHil. There are regions in the AF phase where the magnetization as a function of 
magnetic field has jumps. 

The second considered model is the special case of the spin ladder with exchange integrals depending on one param- 
eter. The exact ground state of this model is non-degenerated singlet with exponentially decaying spin correlations, 
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FIG. 1. Phase diagram of the 'zigzag' model (|l|). The sohd line is the boundary between the ferromagnetic and singlet 
phases. Circles correspond to the special cases of the model. On the dotted line the ground state is a product of singlet pairs. 



and there is an energy gap. 

It will be shown that proposed form of exact wave function can be generalized to higher dimensions and two 
dimensionalal frustrated spin model with exact ground state can be constructed. 

It .is Jen own .that the exact ground state of some ID and 2D quantum spin models can be represented in RVB 
fornxJ'BESLjEj. The RVB function in the Fermi representation consists of homopolar configurations of electron pair. 
We generalize this two-particle function to include ionic states and denote it as "singlet bond" (SB) function. It 
is natural to try to find electronic models with exact ground state formed by SB functions in the same manner as 
for known spin models. In the paper some models of interacting electrons are presented. The Hamiltonians of these 
models include the correlated hopping of electrons and spin interactions. One of the ID models describes the transition 
point between the phases with and without an off-diagonal long-range order. 

The paper is organized as follows. In Sec. 2 we consider the frustrated spin chain at F-AF transition point and 
describe the exact singlet ground state wave function as well as details of the spin correlation function calculations. 
We discuss the phase diagram of this model and its magnetic properties in the AF phase. In Sec. 3 the special spin 
ladder will be considered. Two-dimensional frustrated spin model with the exact ground state is considered in Sec. 4. 
Sec. 5 is devoted to the construction of the electronic models with SB type of wave function. The results of paper is 
summurized in Sec. 6. 



II. ZIGZAG SPIN MODEL. 
A. Zigzag spin model at F-AF transition point. 

Let us consider s = ^ spin chain with nearest- and next-nearest neighbor interactions given by the Hamiltonian 

M M N ^ 

H=- ^(S2.-l • S2, - -) + J23 • S2.+1 - -) + ./l3 ^(S. • S,+2 - (1) 

i—1 i—1 2—1 

with periodic boundary conditions and even number of spins N — 2M . This model is equivalent to a ladder model 
with diagonal coupling, so called a 'zigzag' spin ladder. 

If J23 < 1, then the ground state of (||) is ferromagnetic (singlet) at 5 < (<? > 0), where S = J23 + i-'i^J^^ (Fig0)- 
The equation S — defines the line of transition points from the ferromagnetic to the singlet state, when energies 
of these states are zero. The model (|^) along this line is given by the Hamiltonian depending on the parameter x 
{x > -1/2): 

M I 2x ^ 1 ^ 1 

H = - ^(82,-1 • S2. - -) + ^(82. • 82.+1 - -) - x^(8. . 8.+2 - -), (2) 

i=\ i—1 i—1 

with periodic boundary conditions. 
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We note that the Hamiltonian (||) has a symmetry: its spectrum coincides with the spectrum of H{x) obtained by 
the following transformation 

H{x) = H{-x --), --<x<0 

^ ' 2x+l ^ T 2 



This transformation permutes the factors at the first and the second terms in the Hamiltonian 



the symmetry it is sufficient to consider the range a; > ^ 



4- 



First, we will show that the ground state energy of is zero. Let us represent the Hamiltonian 
Hamiltonians of cells containing three sites 



Thus, due to 
j) as a sum of 



M 



(3) 



where 



X 1 1 



1 1 X 1 1 

= -2(^22+1 • S2i+2 ^ 4) + 2x+l ' ^^*+^ ^ ^ 2;(S2,; • S2j+2 - ^) 



Eigenvalues of each /i„ are 



X X n X 4:x'^ + 2x + l ^ 

A,=A2=0, A3= ^^^2 >0 



We will present a singlet wave function which is the exact one of each /i„ with zero energy and, therefore, it is the 
exact ground state wave function of (|^). This function has a form 

*o = ^^o*, = {Ob\ gi (Si g'2 <S> 93 <E) . ■ ■ <Si On \0b) (4) 

where 

9^ = b+\^)^ + \i)^. gl ^ [b+b ~ 2x) \^) ^ + b \i) ^ (5) 

Here we introduced one auxiliary Bose-particle 6+ (the Bose operators 6+ and b do not act on spin states |t)j and 
and the Bose vacuum |0b). Therefore, the direct product gi g'2 ... g'j^ is the superposition of all possible 
spin configurations multiplied on the correspondiM Bose operators, like b^bbb^ . . . |tiJ,T • ■ •)• -Po is a projector onto 



the singlet state. This operator can be written asl 



L6 



^0 = ^/ da dpi sin7d7 e*"^'e'^•^'e''3^^ (6) 



^0 Jo 



where S"^*^^^ are components of the total spin operator. 

The form of wave function (^) resembles the MP form, but with an infinite matrix which is represented by Bose 
operators. Thereforeyi we have to pick out the (Oti| . . . |0f,) element of the matrix product instead of the usual trace 
in the MP formalismEl't3, because the trace is undefined in this case. The function 5" contains components with all 
possible values of spin 5 (0 < < N /2) and, in fact, a fraction of the singlet is exponentially small at large N . This 
component is filtered out by the operator Po- 

One can easily check that each cell Hamiltonian /i2i-i and /i2i for i = 1, . . . (M — 1) gives zero when acting on the 
corresponding part in wave function - 521-1 ® 52i ® 92i+i and (/ji ® ff2i+i ® 32i+2 (one should take care of Bose 
commutation relations). Since each hi is a non-negatively defined operator, then ^E* is the exact ground state wave 
function of an open chain: 

Hopcn — ^ (^2i-l + h2i) (7) 
1=1 

As mentioned above, the function ^ contains components of all possible valup|[Of total spin S, and, therefore, the 
ground state of the open chain is multiply degenerate. However it can be provenEJIlil that for the cyclic chain (0) only 
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singlet and ferromagnetic components of vj/ have zero energy. Therefore, for cychc chain (|3|) Vfo is the singlet ground 
state wave function degenerate with the ferromagnetic state. 

In particular case, x — — l/4,jjzhen J12 = J23 = — 1 and J13 = another form of the exact singlet ground state 
wave function has been found inEj. It reads 

vI' = ^[i,j][fc,Z][m,n]..., (8) 

where denotes the singlet pair and the summation is made for any combination of spin sites under the condition 
that i < j, k < l,m < n . . .. 

The following general statements relevant to the Hamiltonian (^) can be proved: 

1) . The ground states of open chain described by (|^) in the sector with fixed total spin S are non-degenerate and 
their energies are zero. 

2) . For cyclic chains the ground state in the S — sector is non-degenerate. The ground state energies for 
< S" < M are non-zero. 

3) . The singlet ground state wave function for open and cyclic chains coincide with each other. 



B. Spin correlations in the ground state. 

For the sake of simplicity we show the calculation of the spin correlation function in the symmetric case x = —1/4, 
when the Hamiltonian (|^) takes the form 

N 11^ 1 

H = - 5](S. • S.+i - -) + - ^(S. • - (9) 

i=l i=l 

Since in this case there is one spin in elementary cell, the singlet ground state wave function ^'o can be written in 
a more simple and symmetric form; 

*o = /'o*, * = (Obi gi® 32 SwiOb), (10) 

where 

ff. = &+ IT>. + Mi>. (11) 

One can check that wave functions ( |l0| ) is the singlet ground state wave function with zero energy of Hamiltonian 
(^. Therefore, the equivalence of wave functions (0) and (|l0|) follows from the non-degeneracy of ground state in the 
5 = sector (though functions in Eqs.(|To|) and^) are different). 

Now we calculate the norm and correlation function of the wave function ^0 (jl^)- The norm of the singlet wave 
function 5*0 is 

(fo|*o) = (*|Po|*) (12) 
It is easy to check that the function 5* has S"^ = 0. Then the projector Pg in Eq.(^ takes the form0 

Po^l f smjd^e''-'^'e"'^\ (13) 
^ Jo 

where z — tan ^, z' = sin ^ cos ^ and S^*-^-* are the operators of the total spin. 
Therefore, the norm takes the form 



(^ol^o) - 2 / sin7d7(0a,0,|n(ff.+ e^'^-'e^''^'^5.)|0a,0 



where 

gt e^^'-e^^'^'g^ = («+ (T.I + « {U\ ) e-^.~e-'^.*(6+ |T), + b |i), ) 
= a^6^ + (1 — z' z)ab + izab^ + iz'a^b, 

and a+ and a are the Bose operators. Thus the norm can be rewritten as 
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(*o|vl/o) = i rsm7d7(0|G^|0), 



where |0) = |0a, Oh) is the Bose vacuum of a"*" and 6+ particles and 

G = u{a+h+ + ah) + iv{ab+ + a+6), 

where u = cos ^ , u = sin ^ . 



Let us introduce the auxihary function P{C)'- 



PiO = (0|e««|0) 



then 



|0) = 



N 



The function P(^) can be easily foundEl: 



Integrating Eq.(|l^) over 7, we obtain 



\J\- V? sin^ ^ 



(*o|*o) = T, I sin7d7 



Thus, finally, we arrive at 



(*o|*o) = 2^^ (^tani 



=0 rf^^ Vc««'(i) 

4(2^+2-1) 



iV + 2 



P 



N+2\ 



Here Pat are the Bernoulli numbers. 

To calculate the spin correlators we need to introduce operators 



G, = g+ e'''^^ e'^'^^hs^'g., = u{a+b+ - ab) + iv{ab+ ~ a+b), 



G+ = gf e*"^. e*" *V St g^ = ua+h + ivab, 



G- = 5+ e''^^ e"- ^- S^gi = uab+ + iva+b+ 
Then, the correlator (Si • S;+i) will be defined by 



(M/oISi -S^+il^o) ^ ^ J\injdj {0\ \g,G'G,G 



Iri _^ _Q qIq G^~ 



(since (0|G |0) = 0). 

The expectation values in Eq.dlS) can be represented as 



)| G.G'GzG^"'"^ |0) = 
G+G'G_G^-'-2|0) = 



9' d 



Ql qN-1-2 



After a procedure similar to that for the norm and the integration over 7, we obtain 



I aN-l-2 



(*o|Si -Si+il^-o) = 



a' d 



I f)rN-i-2 



3 cos(e - C) 
8cos4(«±«) 



<=0 



It can be shown that in the thermodynamic limit, Eqs.(|l8D and (|2l|) result in 




FIG. 2. The thick sohd hne is the boundary between the ferromagnetic and singlet phases. Circles correspond to the special 
cases of the model. The dotted line denotes the heuristic boundary between gapped and gapless phases. The dashed line is the 
boundary of the region with multimagnon bound state. 

(S, • S,+i) = i COS (^^^ (22) 

So, we reproduce the result obtained inE10 that in the thermodynamic limit a giant spiral spin structure is realized, 
with the period of the spiral equal to the system size. 

For the general case of model (||) the calculation of the singlet gmpjd state correlation functions can be performed 
by the similar way. The final result in the thermodynamic limit isEJO: 

(S„S„+2i) = ^ cos 

(S„S„+2;-i> = \ cos (^^^^^^ + (-1)"A^) (23) 

The latter equations mean that the long-range double-spiral order exists in the singlet ground state of Hamiltonian 
(^. The pitch angle of each spiral is ^ and there is a small shift angle A(/3 = 2;r(4a+i) ]2g^^.^ygg-^ them. This shift 
angle reflects the fact that the unit cell contains two sites unless x — —\. 

We note that for special values of model parameter x = ^,1,|,2... Eqs.(|2|) are not valid and spin correlations 
decay exponentially. These cases will be considered in Sec. Ill (see Eqs.(^)). 




C. Phase diagram of 'zigzag' model 

So far we considered the model (^ at the transition line from ferromagnetic to antiferromagnetic state. Now we 
discuss the phase diagram of this model. The exact ground state in the AF phase is generally unknown. But it is 
interesting to note that the ground state on the line J13 = —1/2 is the product of singlets on ladder diagonals (2,3), 
(4,5), ... as in the point J13 = — J23 = —1/2 on the transition line. The spectrum of (|^) on the transition line is 
gapless. There are some regions on the plane (J13, J23) which were studied by different approximation methods. 

First let us consider behaviour of the system in the transition region near the 'symmetric' point x = —j, where 

Jl2 = J23 = -l, Jl3 = l + 6 (24) 
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It would be mentioned that several copper oxide compounds are described by this model. These compounds contain 
CuO chains and Cu — O — Cu angle 9 is near 90°Ll. In this case the usual antiferromagnetic super-exchange of two 
nearest-neighbor magnetic Cu ions is suppressed and the exchange integral J12 is negative. On the other hand the 
next-nearest-neighbor interaction J13 between Cu ions does not depend on 9 and is antiferromametic. The estimation 
of the ratio J13/J12 for LaQCa^Cu2AOii and Li2Cu02 gives S = 0.11 and S — 0.37 respectivelylD. 

Though the model (|2^ ) is not exactly solvable in this case, its properties for (5^1 can be studied. The consideration 
is based on the classical approximation. In this approximation the ground state spin structure is the spiral with the 
period ^ S~^^^ and the ground state-ejuergy E — ~2N6^. Using this approach the regular expansion in powers of the 
small parameter S can be developedllalU. The second-order quantum corrections coincides with the classical energy. 
The calculation of higher orders of the perturbation theory in S leads to infrared-divergent integrals and it is necessary 
to sum them in all orders to obtain the contributions proportional to 5^/^, 6^ etc. In particular, the ground state 
energy calculated up to terms ~ S^^^ is 

Ea = -mS^ + i.UNS^/"^ 

The excitation spectrum is gapless and has a sound like behaviourEl. Corresponding sound pvelocity is w = iS^^^. 

The study of the dependence of the ground state energy E{m) on magnetization m — 5z/iVE3 at (5 ^ 1 shows that 
< in a finite interval of values of m. This implies the thermodynamic instability of the uniform state against 
the phase separation. This instability arises due to existence of multimagnon bound states. Energy of n-magnon 
bound state (n ^ 1) is e„ = — nfe, b — |£i| -I- ef,, where ei is the energy of one-magnon state and is binding energy 
per magnon. For the model (24) ex = —8(5^ and tb ~ J^/^. As a result the function E{m) has the forrrE3: 

E{rn) = E'o + 27r(5^/^TO^, m < rric ^ —VS 

i?(m) = -^(1 - 2m), mc<m<^ (25) 

At m > rric the system is in a two-phase state consisting of the ferromagnetic (to = 1/2) phase and the phase with 
TO = TOc. According to (p5| ) the magnetization as a function of magnetic field h has a jump frorp-p, = TOc to to = 1/2 
(metamagnetic transition) a.t h = he = b. These values TOc and he are close to those obtained inEZI by extrapolation of 
finite cluster calculations. We stress that the jump of the magnetization exists if there are multimagnon bound states 
leading to the linear dependence of E{m) for some region of m. ___ 

As for properties of the model far from the transition point (5 = 0, now-days there is no cleat-understanding3't2lE2l, 
but we believe that the spectrum remains gapless. For the case i5 ^ 1 this fact is predicted inEil. As for the jump in 
m(h) it is unknown if it remains at 5 ~ 1. Numerical calculation£3 show that, at least, there is an abrupt increase 
of magnetization at some critical value of magnetic field. 

This perturbation theory can be generalized for the vicinity of any point on the transition line. It gives the similar 
properties of the system in the region J23 < 0, but it diverges at J23 > 0. |— ■ 

At Ji3 = and J23 > the model (|l]) reduces to the alternating Heisenberg chain studied ii£3. The lowest 
excitation is the triplet and there is the gap. At J23 = and J13 > the model (|l|) reduces to the spin ladder with 
antiferromagnetic interactions along legs and the ferromagnetic interactions on rungs. It is evident thatnthere is a gap 
at Ji3 <C 1 (in this case the model is equivalent to the spin S = 1 Heisenberg chain). It was shown inE3 that the gap 
exists at J13 ^ 1. nn 

The region J23 > 0, J13 < was studied by different methods inESEa. The exact diagonalization of finite chains 
shows the gap A in the excitation spectrum which is closed on the transition line. The presence of singlet-triplet 
gap is also confirmed by calculations with use of variational wave function in MP form ( |39| ) whcii-alLmatrix elements 
are considered as variational parameters. The MP variational function gives very good accuracyu'ES, besides it gives 
exact ground state results for special points on the transition line and for the line J13 = — ^ (dimer state) . 

Another important property of the model with J23 > 0, J13 < is the existence of the multimagnon bound states 
for definite region of the values J23 and J13. This region is shown on Fig.l. On the boundary of this region the 
multimagnon bound states disappear. As it was noted before the existence of the bound states leads to the linear 
regions on the dependence E{m) and to the metamagnetic transition. The dependence E{m) and m{K) was found 
using MP variational function. The calculations show that out of the bound state region the function m{K) is typical 
for the gapped antiferromagnet, i.e. it is monotonically increasing function (Fig.^). In other words, m{h) = for 
/i < A and m{h) — 1/2 for h > \ei\. In the bound state region the magnetization curve has the jump. It is interesting 
that inside the stripe —\< J13 < ~^ this jump takes place immediately from zero to the maximal value to = ^ 
(Fig.^) at he = |e(0)| (e(0) is a ground state energy per site). Such behaviour of m{h) is explained by the fact that 
inside this stripe A > |£(0)| and the transition from the singlet state to the ferromagnetic one comes to pass through 
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m 



h 



FIG. 3. Dependence of magnetization m on magnetic field h for the case: (a) outside bound state region; (b) in the bound 
state region, but outside the stripe — ^ < Jii < (c) inside the stripe —\< Jis < 
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FIG. 4. The two-leg spin ladder. 



states with intermediate values of spin. In the bound state region outside the stripe the magnetization curve has a 
form as shown on Fig.(|). In this case the critical field he — A. 

Thus, summarizing al above, we expect that the phase diagram of the model (^) has the form shown in Fig.(H). 



III. SPIN LADDER MODEL 



Let us consider more general spin ladder model, for which 'zigzag' model is a particular case. So, we consider 
the cyclic ladder model containing N = 2M spins s — ^ (Fig.^. The proposed form of wave function can be 
generalized for a ladder model as follows: 



*o = -Po*, * = (0| 51 (8)^2 «)...«) gA/|0), 
where each gi corresponds to the ith rung of the ladder: 

9^ = -b+{2x - b+b) ITT), + lii), + {b+b - x) \n + it), + y\n- it), 

where x and y are two parameters of the model. It is obvious that the ladder wave function ^ 
original function at y = x. 

Now we will construct the Hamiltonian for which 5'iaddcr is the exact ground state wave function. This Hamiltonian 
describes two-leg s = 5 ladder with periodic boundary conditions (FigH) and can be represented in a form 



(26) 
(27) 

|27D reduces to the 



H 



N 
i=l 



(28) 



where describes the interaction between neighboring rungs. The spin space of two neighboring rungs consists of 

six multiplets: two singlet, three triplet and one quintet. At the same time, one can check that the product gi (g) gi+i 
contains only three of the six multiplets of each pair of neighboring rungs: one singlet, one triplet and one quintet. 
The specific form of the singlet and triplet components present in the product gi (g) gi+i depends on parameters x 
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and y. The Hamiltonian can be written as the sum of the projectors onto the three missing multiplets with 

arbitrary positive coefhcients Ai, A2, A3: 



K,+,=Y^XuPi:'+\ (29) 
fe=l 

where P^'*^^ is the projector onto the missing multiplets in the corresponding cell Hamiltonian. 

The wave function ( p6|) is an exact wave function of the ground state of the Hamiltonian with zero energy, 

because 

=0, z = l,...iV-l (30) 

and Ai, A2, A3 are the excitation energies of the corresponding multiplets. 

So, \1/ is the exact ground state wave function with zero energy for the total Hamiltonian of an open ladder 

M-l 

Hopcn — ^ hi^i+i (31) 



ifopcn|*>-0 (32) 

Since the function contains components with all possible values of total spin (0 < < M), then the ground 
state of open ladder is multiple degenerate. But for a cyclic ladder ( |2^ ) only singlet and ferromagnetic components 
of 5" have zero energy. Therefore, for a cyclic ladder is a singlet ground state wave function degenerated with 
ferromagnetic state. Besides, all the statements given for 'zigzag' model (||) are valid for ladder model (p8|). 

Since the specific form of the existing and missing multiplets in the wave function (^6|) on each two nearest neighbor 
spin pairs depends on the parameters x and y, the projectors in ( p9| ) also depend on x and y. Each projector can be 
written in the form 

P^'^ = "^12'' (^1 ' ^2 + S3 • S4) + j}3''(Si • S3 + S2 • S4) + ■ S4 + J23''^2 • S3 

+ jf )(Si • S2)(S3 • S4) + ^(Si • S3)(S2 • S4) + (Si • S4)(S2 • S3) + (33) 

and this representation is unique for a fixed value of the parameters x and y. 

Substituting the above expressions for the projectors into Eq. (|29|), we obtain the general form of the Hamiltonians 
hi^i^i. Inasmuch as the Hamiltonians have the same form for any i, it suffices here to give the expression 

for hi^2'- 

hia = Ji2{Ai2 + ^34) + Ji3(^i3 + ^24) + JuAii + J23A23 

+ J1A12A34 + J2A13A24 + J3AUA23 (34) 



where 



Ai-i — Si • S7 



1 



and all exchange integrals depend on the model parameters and the spectrum of excited states Ji — Ji{x, y, Ai, A2, A3) 
as follows: 

^ _ A2 A3V_i 7 _ ^2 A3 (2?; -1)2 

^12 ~ TT ' TT ^ 9 I 1 ' "^23 



2 2 iy^ + 1' ' 2 2 4y2 + 1 ' 

J A2 A3 - 1 A2 A3 (2y + 1)2 

"~ 2 2 4y2 + l' -^14- 2 2 4y2 + l ' 

TOT X y^-x^{x + lf T 07 . o\ 2/^ + y^a:(a; + 1) 



-3y4 + a;2(a; + l)2' " ' 3?/4 + a;2(a; + 1)^ 

y** - y^x^x + 1) 
-3y4 + a;2(a; + l)2 



7 7 . 7 . o\ v^-y^x{x^\) 

^3 = -/m + J23 + 2Ai— — — — — — — - (35) 
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(one should keep in mind that only positive \i can be substituted to these expressions). The model (^5|) has evident 
symmetry: the change of sign of y is equivalent to renumbering of sites 1^2, 3 <-!■ 4 . . .. Therefore, we will consider 
only the case y > 0. 

In general, the Hamiltonian contains all the terms presented in (|3^), but we can simplify it by setting, for 

example, J2 = J3 — and solving equations ( p5| ) for Ai, A2, A3. All Afc turn out to be positive in this case for any x 
and y except two lines y = and x = —1/2, where ground state is multiple degenerated. The Hamiltonian in 
this case takes the form 

hl^2 = Jl2(Ai2 + A34) + Jl3(^13 + ^24) + Jl4^14 + ^23-423 + ^1^12^34 (36) 

Ji2 = -x{x + 1) + - ^, Ji4 = x{x + 1) - - y, 

Ji3 = -x{x + l)-y\ J23=x{x + l)-y^ + y, = M^±i)_^ _ 1 

y2 

The calculation of the norm of ( p6| ) and the singlet ground state correlation functions can be performed in the 
similar way to the corresponding calculations for the case y = x = —1/4. Therefore, we give here the final result for 
spin correlation functions at iV — > oo 

(S„ • S„+2/) = ^ cos , 

1 /47]-/ \ 

(S2„-i • S2„+2i) = - cos ( — - A</7 j (37) 

These equations mean that the spiral on each leg with pitch angle ^ is formed and the shift angle between spirals 
on the upper and the lower legs is A<p = At y = 0, when spins on each rung form a local triplet, the shift angle 
vanishes and the spirals on both legs become coherent (we note that the shift angles in Eqs.(23) and (^7|) are defined 
in a different way). 

Thus, there is just one full rotation of the spin over the length of the ladder, independent of the size of the system 
and for fixed ^ ^ iV at — > oo two spins on the ladder are parallel. 

We emphasize that the spin correlation function (S^ • Sj) does not depend on the choice of Afe for a fixed parameters 
x, y, because the ground state wave function of the three-parameter set of Hamiltonians (Isl) is the same. 



A. Special cases 

There are special values of the parameter x for which Eqs.(|3^) are not valid. For cases of integer or half-integer 
X = j, which correspond to the specialjiases of the model (^5[), in Eq.(p7|) one can easily recognize the Maleev's boson 
representation of spin S — j operatorsEj: 

S+ = b+{2j -b+b), S-=b, S' = b+b-j (38) 

Generally, the wave functions (Q) and (^) resemble MP form but with an infinite matrices represented by the Bose 
operators. However, in accordance to Maleev's representation in the special cases the infinite matrices formed by the 
Bose operators 6+ and b can be broken off to the size n = 2j + 1 and wave function (|2^) is reduced to the usual MP 
form 

*o = (gi ® 32 ® ■ • • <^ 9n) , (39) 

where g = —xT + yS is the n x n matrix describing states of spin pair on corresponding rung of the ladder. Singlet 
state matrix is 

S^I\s) (40) 

where / is unit matrix and \s) is the singlet state. Triplet state matrix T is expressed by Clebsch-Gordan coefficients 
Cmi,m2 = ((l,"ii) (j,m2) II (i,mi -f- 1712)) as follows: 
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FIG. 5. Stripe spin structure on the ladder model. 
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Co,,-i|0) . 
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Ci,-, |i) 



(41) 



C_i,_,+i 1-1) Co,-j |0) / 



where \a) is the triplet state with = a. 

Exact calculation of the correlators in the thermodynamic limit is performed using standard transfer matrix tech- 
nique and results in 



(S1S2) 

(S,;S,i+2/) 
(S2i+lS2i+2/+2) 
(821+2821+2/+!) 



x{x + 1) - 

x{x + 1)(4?;2 - 1) fuj-l 
x{x + l){2y-lf fuj-l 



i~i 



4w2 



OJ 



x{x + l){2y + If fuj- 1 



where 



4w2 



u) = x{x + 1) + 



OJ 



(42) 



In the particular case of zero singlet jiueight, y = 0, when spins on each rung form a local triplet, correlation functions 
( p^ ) coincide with those obtained incl'Ej. 
According to Eqs.(^) the spin correlations have an exponential decay and the correlation length Tc is 



= 21n"i 



U! - 



(43) 



In particular, for special points of 'zigzag' model (||) with J23 — 2x+i "^13 ~ ~^ correlation length = 

-21n-i ( 1 - -pT^). 

y x{2x+i) J 

The correlation length rcdiverges when x ^ 00 or y ^ 00. In these cases the singlet ground state has collinear 
or stripe spin structure, i.e. spin-spin correlations are ferromagnetic along legs and antiferromagnetic between them 
(Fig.^), with a magnetic order m: 



(S,;S 



1^4 + 2// 



(S,S 



xy 



x^ -f y2 



When y = X (the 'zigzag' model) the magnetic order is equal to the classical value m = 1/2. 

We note that the wave function ( p6|) show double-spiral ordering for all values of x and y excluding the special lines 
X = j. The crossover between spiral and stripe states occurs in the exponentially small (at N 00) vicinity of the 
special lines. 
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B. Spectrum of the model 



Generally, the excitation spectrum of model (^^6|) can not be calculated exactly. It is clear that this spectrum is 
gapless because, for example, the one-magnon energy is ^ N^^ dX N ^ co. Moreover, the lowest singlet excitation 



is gapless as well. For model (^|36) lying on the special lines this fact can be established from the following 
consideration. As mentioned above, the crossover between spiral and stripe states occurs in the exponentially small 
vicinity of the special lines. It means that two wave functions ^'o(j) and + 6) corresponding to points {x = j, y) 
and {x = j + S, y) respectively are almost orthogonal at (5 ~ . Hence, we can consider a spiral wave function 
*o(j + S) as a variational function of the excited singlet state at special point x — j. The energy of this excited state 
at special point x ~ j is 



*oO' + -5)1 H{j + 5)- 5 \^o{3 + 5)J - const. • S 

Thus, on the special lines the stripe ground state of considered model is asymptotically degenerated with an excited 
spiral singlet state at the thermodynamic limit. It is not clear if the degeneracy is exponentially large or not. This 
consideration is valid for any integes or half-integer x — j, but it is not valid for the parameters x ^ j, which are 
out of special lines. We performedcJ the numerical diagonalization of finite ladders for various parameters x and y 
and found that the exponential degeneracy possibly takes place for all parameters x and y, but we can not confirm it 
strictly. |— .|_. 

It is interesting to note that singlet wave function (|6|) can be also represented in a special recurrent foriroO 

*o = ^o*M, (44) 



*M = {st + Vist + V2S^ ...+ !^2S^)(S3 + i^lS4 . . . + ^28%) ■ ■ ■ 

{4n-l + ^l4n ■■■ + ^^24) ■ ■ • (4-1 + ^^14) lii • • • i) (45) 

where s'l is the s — ^ raising operator. Eq.(^) contains M operator multipliers and the vacuum state \[l ... [) is 
the state with all spins pointing down. The function 'I'm is the eigenfunction of with 5^ = but it is not the 
eigenfunction of S^. Pq is a projector onto the singlet state. Two parameters vi and V2 in wave function (^) is 
connected with parameters x and y in ( |2^ ) as 

1 + X — y 1 
1 + X + y \ + X + y 

The normpafj-the wave function (^) and expectation values can be also calculated with use of recurrent technique 
developed iniJO. Certainly, it gives the same expressions (K^E2h for spin correlation function. 



C. Antiferromagnetic ladder model 

Now we should consider in particular the special case x — 1/20. In this case the product gi (E) gi+i contains only 
one singlet paiid one triplet and does not contain quintet. Therefore, in this case the cell Hamiltonian can be written 
in the formO 

4 

H = J2 h^+i^ K^+i = E ^kPt:'^\ (46) 

i k=l 

where Pl'^~^^ - is a projector onto the quintet state. If all \k > 0, the ferromagnetic state has positive energy E = MA4 
and wave function 'I'o ( p6|) is now non-degenerate singlet ground state wave function for the Hamiltonian (46), while 
for Hamiltonian ( p^ , p9| ) 'I'o is also exact ground state but degenerate with ferromagnetic state. Using a freedom in 
choice of we can exclude all four-spin interactions in 

hi,i+l — J12 (S2i-1 • S2i + S2i+1 • 82^+2) + Jl3 (822-1 • 821+1 + 82i • 821+2) 

+ Jl482i-1 • 82^+2 + J23^2i ■ 82^+1 + C (47) 
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A(y) 




FIG. 6. Singlet-triplet gap of the model (g8|) as a function of the parameter y. The circles denote the results of the 
extrapolation of an exact finite-chain calculations. The solid line represents the dependence A(j/) given by Eq.(^o|). 

and all exchange integrals Jy depend on one model parameter y < 3/2 (this inequality is necessary to satisfy > 0) 



Ji2 = 2 (V - 1), Ji4 = -y(3 + 2j/)(2y - 1)2, 

Ji3 = ~2v\Ay^ - 1), J23 = 2/(3 - 2y)(2y + 1)\ 



^ = + 7 



(48) 



It follows from Eg. ([42|) that the ground state has ultrashort-range correlations with rc ~ 1. For example, rc{y — 
0) = 21og^^ 3, which coincides with correlation length of AKLT model. But dX y — ^ all correlations are zero except 
(S2i • S2i+i) — — J- It implies that at y = ^ model ( ^6[]4^ ) has a dimer ground state. 

The value (w— 1) changes the sign aty—^ and as it follows from Eqs.(^^ the correlators show the antiferromagnetic 
structure of the ground state at ^ < y < | while at < y < i there are ferromagnetic correlations inside pairs 
(1, 2), (3, 4), . . . and the antiferromagnetic correlations between the pairs. 

The Hamiltonian (^ -48) of the cyclic ladder has a singlet-triplet gap A for finite N. It is evident that for y = 5 



the gap exists for iV — > cx) and A(i) = 4. The existence of the finite gap at the thermodynamic limit in the range 
< y < I follows from the continuity of the function A(y). It is also clear that A(y) at iV ^ cxd vanishes at the 
boundary points y = and y = § when the ground state is degenerate and there are low-lying spin wave excitations. 

Unfortunately, the method of the exact calculation of A(y) in the thermodynamic limit is unknown. For the 
approximate calculation A(y) we use the trial function of the triplet state in the form 



E 



+ ikn 

sle 



(49) 



The trial function '^t gives A(y) which at ^ 00 has minima at fc = tt and A: = OforO<y<i and ^ < y < |, 
respectively 



A(y) = 
A(y) = 



32y2(4y2 



1) 



4y2 3 ' 
128y2 
(4y2 + l)(4y2 + 3)' 



0<y<l 



1 3 
- < y < - 
2^2 



(50) 



The dependence of A(y) given by Eq.(|50|) is shown on Fig.^ together with the results of extrapolations of exact 
finite-chain calculations. Both dependences agree very well for y < 5. However, A(y) given by Eq.(50) is not zero at 



y = |, while numerical calculations fit the dependence A(y) ^ y | — y at y ^ |. 

We note that the trial function of the type ( p9[ ) gives the value 0.7407 for Ihe singlet-triplet gap in the AKLT model. 
This estimate is close to the value 0.7143 obtained by another approach inEJ. 

For other special cases one can also construct Hamiltonians for which V&o is non-degenerate singlet ground state 
wave function. But in these cases one have to introduce more distant interactions. 
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IV. VALENCE-BOND-STATE MODELS 



A. One-dimensional model 



We studied previously a one-parameter ladder model (|46|-[18D with non-degenerate singlet ground state. The exact 
ground state wave function of the cyclic ladder was written in the MP form Now we write the wave function ^'o 
in a form more suitable for subsequent generalization to other types of latticeH. 

We consider a ladder of TV = 2M spins 1/2. The wave function of this system is described by the iVth-rank spinor 

^ ^^\t^t....T^ (51) 

where the indices A, /i, i', . . . , t = 1, 2 correspond to different projections of the spin 1/2. 

We partition the system into pairs of spins located on rungs of the ladder. The wave function can then be written 
as the product of M second-rank spinors 

^^<jf^M(i)^!^P(2)...*'"^(M). (52) 

We now form a scalar from Eq. (^2|), simplifying the latter with respect to index pairs: 

-1'\(1)*^(2)...^''^a(M). (53) 

Here subscripts correspond to the covariant components of the spinor, which are related to the contravariant compo- 
nents (superscripts) through the metric spinor 

3A, = .9'^=(_°1 ly (54) 

The scalar function (^3|) can thus be written in the form 

- vl/V(i)g^^^vl/-P(2)5^, . . . ^''^{M)grx. (55) 

The scalar function 'i's obviously describes to the singlet state. 

The second-rank spinor describing the pair of spins 1 /2 can be written in the form 

=ct*f''+c,*^'^, (56) 

where '^/^^ and '^^^ are symmetric and antisymmetric second-rank spinors, respectively, and ct and Cs are arbitrary 
constants. We know that the symmetric second-rank spinor describes a system with spin S* = 1, so that the pair of 
spins 1/2 in this case forms a triplet. If is an antisymmetric second-rank spinor reducible to a scalar multiplied 
by g\^l, the spin pair exists in the singlet state. Consequently, the ratio of the constants c* and Cs determines the 
relative weights of the triplet and singlet components on the pair of spins s = 1/2 and is a parameter of the model. 
In particular, for Cg — the wave function ( p6| ) contains only a triplet component, and for ct — it contains only a 
singlet component. 

We note that the wave function ( ^3| ) has MP form (|3^) with the matrices gi representing a mixed second-rank 
tensor: 



^^Ti + iT), lii), \ 1 



2 1 

Y -ITT). -2'^^ + ^^^ 



- ^csiin + iT). /, (57) 



where / is the unit matrix. 

We now choose a Hamiltonian H for which the wave function ( |55| ) is an exact ground-state wave function. To do 
so, we consider the part of the system (cell) consisting of two nearest neighbor spin pairs. In the wave function ( [ss] ) 
the factor corresponding to the two spin pairs is a second-rank spinor: 

vl/^^(*)g^,vl/-P(* + l). (58) 

In the general case, therefore, only two of the six multiplets forming two pairs of spin 1/2 — one singlet and one 



triplet — are present in the function (58). Inasmuch as four spins 1/2 form two singlets and three triplets, the specific 
form of the singlet and triplet components present in the wave function (^8|) depends on the ratio Cg/ct- The cell 
Hamiltonian acting in the spin space of nearest neighbor spin pairs can be written as the sum of the projectors onto 
the four missing multiplets with arbitrary positive coefficients (^6|). As mentioned above, the general form of the 
Hamiltonian (Hq) can be reduced to a more simple form (|47|) with Cs/ct = 2y. 



Thus, the singlet ground state wave function of the model (Mq) can be also written in a spinor form (55) 
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FIG. 7. Two-dimensional lattice on which the SDin model is defined. 




FIG. 8. Graphical correspondence of the model wave function. The indices of the site spinors depend on the site index (not 
shown in the figure). 



B. Two-dimensional model 



Now we consider an iV x TV-site square lattice with cyclic boundary conditions. We replace each site of the lattice 
by a square (Fig.0) with spins s = 1/2 at its corners, making the total number of spins equal to AN'^. To avoid 
misunderstanding, however, from now on we continue to refer to these squares as sites. The wave function of the 
system is described by the product of fourth-rank spinors 



(n). 



(59) 



By analogy with ( |55| ) , from Eq. (|59[) wc form the scalar 



(60) 



where a and b are unit vectors in the x and y directions. 

The singlet wave function ( |60| ) is conveniently identified graphically with a square lattice, each site corresponding 
to a fourth-rank spinor ^^■^^^'^P (whose form is identical for all sites), and each segment linking sites corresponds to a 
metric spinor .g^^ (Fig.|). 

To completely define the wave function (|6^), it is necessary to know the form of the site spinor 

■qjXf.i'P rpj^g specific 

form of the fourth-rank spinxu: ^i-^'^'^p [and, hence, the wave function (|60|)] describing the system of four spins s = 1/2 
is governed by 14 quantitiesBil, which are parameters of the model. 
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FIG. 9. Lattice sites associated with interactions -ffi,2 and -ffi,3. 

We now choose a Hamiltonian H for which the wave function (|60| ) is an exact ground state wave function. As in 
the one-dimensional case, we seek the required Hamiltonian in the form of a sum of cell Hamiltonians acting in the 
space of two nearest neighbor spin quartets: 

H = ^ Hn,n+a + ^ Hn,n+h- (61) 
n n 

The first term in Eq. (|^) is the sum of the cell Hamiltonians in the horizontal direction, and the second term is 
the same for the vertical. The cell Hamiltonians along each direction have the same form, but the "horizontal" and 
"vertical" Hamiltonians differ in general. In the following discussion, therefore, we consider only the Hamiltonians 
Hi^2 and i/i,3 (Fig.^), which describe interactions of "sites" in the x and y directions, respectively. 

For the wave function ( |60[ ) to be an exact eigenfunction of the Hamiltonian H, it is sufficient that the sixth-rank 
spinors 

^A1M1..1P1^)^A3M3^3P3(3)^^^^^^ (62) 

be eigenfunctions of the corresponding cell Hamiltonians Hi 2 and iJi,3. 

In general, when the site spinor vp-^^f^p j^gt symmetric with respect to any indices, the possible states of two 
quartets of spins s = 1/2 consist of 70 multiplets. A wave function represented by a sixth-rank spinor contains only 
20 of them. Accordingly, the cell Hamiltonians Hi 2 and Hi 3 can be represented by the sum of projectors onto the 
50 missing multiplets: 

50 50 

Hi.2 = ^ \kPt\ i?i,3 = ^ ^^kPt^ (63) 

fe=l k=l 

where the positive constants Afc and fj,k are the excitation energies of Hi,2 and i/i^s, and the specific form of the 
projectors depends on 14 model parameters. 
Inasmuch as 

ifn,n+a|*s) = 0, i?„,„+b|*s) = 0, (64) 

for the total Hamiltonian (|6|) we have the expression 

H\^s)^0- (65) 

Consequently, 'i!s is the ground-state wave function of thfi-tptal Hamiltonian H, because it is a sum of nonnegative 
definite cell Hamiltonians. Also, it can be rigorously provedEll that the ground state of H is nondegenerate. 

As mentioned above, the specific form of the projectors depends on 14 model parameters, and in general the cell 
Hamiltonians (|6^), expressed in terms of scalar products of the type Si ■ Sj, (si • Sj){sk • s;), etc., have an extremely 
cumbersome form. We therefore consider a few special cases. 
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FIG. 10. Pattern of independent singlet pairs (double lines). 



When the site spinor vp-^Af^p jg ^ symmetric fourth-rank spinor Q^^^'^p (corresponding to the two-dimensional AKLT 
modellj) , only the quintet component out of the six multiplets on each spin quartet is present in the wave function 
(|60|). The sixth-rank spinors ( |62| ) are symmetric with respect to two triplets of indices and, hence, contain four 
multiplets with 5 = 0, 1, 2, 3 formed from two quintets. Consequently, the cell Hamiltonian (-ffi,2 and Hi ^ coincide 
in this case) has the form 



66 



fc=i 



If we set Xk — 1 (k = 1, 66), we can write Eq. ( p6[) in the form 

i?l,2 =^4(81+82)+ [l-P2(Sl)F2(S2)], 



(66) 



(67) 



where Si is the total spin of the quartet of spins s = 1/2 on the ith site. Si = Si{i) + S2{i) + s^ii) + §4(1), and P;(S) 
is the projector onto the state with spin S = I. 

If the four spins s = 1/2 at each site are replaced by a single spin S — 2 and if the wave function ( |60| ) is treated 
as a wave function describing a system of N'^ spins S" = 2, the second term in the Hamiltonian (67) vanishes, and we 
arrive at the Hamiltonian of the two-dimensional AKLT model: 



i/,2 = MS, + S2) = -82 + ^^(8,. S2)^ + ^(S, . S2)3 + 



(81-82)^ 



(68) 



Another interesting special case is encountered when the system decomposes into independent one-dimensional 
chains. This happens if the site spinor ^'^'^'^p reduces to a product of two second-rank spinors, each describing two 
spins 1/2. For example, 



*^^''''(S1, S2, S3, Si) = ^^"{S,, S3)ip^Pis2, S4). 



(69) 



In this case the Hamiltonians i?i,2 and i?i.3 contain interactions of four rather than eight spins 1/2 and have the 
form (11). 

The simplest case is when the site spinor vj/-^'^"'' is a product of four first-rank spinors: 



^^f^-Pisi, S2, S3, S4) = 'p\si)ip^{s2)ip''{s3)ip''{si). 



(70) 



Now the system decomposes into independent singlet pairs (Fig.|l^), and the total Hamiltonian of the system has the 
form 



2J 



(71) 



where s.^ and Sj are the spins forming the singlet pairs. 
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C. Spin correlation functions in the ground state 



We now look at the problem of calculating the norm and the correlation function of the model described by the 
wave function (|60|). The expression for the norm of the wave function G = s) has the form 

n 

= LK"":':::"::: - n^""/5"""+»/5„,„ a„ a = {i, 2, 3, 4}, (72) 



where Rapap isa4x4x4x4 matrix. 

According to the selection rules for the projection of the total spin 5*^, only 70 of the 256 elements in the expression 

M 1^ P[n) '^''^P'^P {n)^ are non-vanishing. Consequently, the matrix R also contains at most 70 elements. If we 

regard the elements of R as Boltzmann vertex weights, the problem of calculating the norm reduces to the classical 
70-vertex model. Since the exact solution for the 70-vertex model is unknown, numerical methods must be used to 
calculate the norm and the expectation values. 

To calculate the above-indicated expected values, we carry out Monte Carlo calculations on 20 x 20-site lattices. 
As mentioned, the ground-state wave function of the model depends on 14 parameters and, of course, cannot possibly 
be analyzed completely. We confine the numerical calculations to the case in which the spinor 

^\pvp depends on one 

parameter a: 

■^xpup = cos a • Q^^'^P sin a • (^A^p'^p - Q^p'^p^ , (73) 
where a G [ — 7r/2; 7r/2] , the spinor Q^p^p is symmetric with respect to all indices, and 

A^P'^P = ^^{si)^P{s2)^^s^)^P{si). (74) 

In this case we have a one-parameter model with two well-known limiting cases. One corresponds to a = 7r/4, for 
which ^^P'^P = A^P'^P ^ and the system decomposes into independent singlet pairs (Fig. |l^); the other limiting case 
corresponds to a = (our model reduces to the two-dimensional AKLT model in this case, the spins at each site 
forming a quintet). 

In the given model there are four spins s = 1/2 at each site, and the enumeration of each spin is determined by the 
order number of the lattice site to which it belongs and by its own number at this site. The spin correlation function 
therefore has the form 

/,,(r)-(s,(n).s,(n + r)). (75) 

In determining the spin structure of the ground state, however, it is more practical to consider the more straight- 
forward quantity F{y): 

4 

F(r)= ^ (s,(n).s,(n + r)) = (S(n).S(n + r)). (76) 

ij = l 

The function F{v) is left unchanged by a change of sign of a. We note, however, that only the total correlation 
function, and not fijir), possesses symmetry under a change of sign of a. This assertion is evident, for example, in 
Fig. [l|, which shows the dependence of /31(a) on a as an illustration. 

The correlation function decays exponentially as r increases, differing from the one-dimensional model in that the 
pre-exponential factor also depends on r . Figure |lj shows the dependence of the correlation length Vc on the parameter 
a. The correlation length is a maximum at the point a = (two-dimensional AKLT model), decreases as a increases, 
and at a = 7r/4, when the system decomposes into independent singlet pairs (Fig. |l0| ), it is equal to zero. With a 
further increase in a the correlation length increases and attains a second maximum at a — 7r/2. Like the correlation 
function F{r), the function rda) is symmetric with respect to a. It is evident from Fig. |l^ that the parameter a has 
two ranges corresponding to states with different symmetries. In the range \a\ < n/A the correlation function F{r) 
exhibits antiferromagnetic behavior: 

Fir) cx (-l)'^=^+^»e"l'"l/''% (77) 
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FIG. 11. Dependence of the spin correlation function (s3(l)si(2)) on the parameter a. 




FIG. 12. Dependence of the correlation length on the parameter a. 



whereas the spins at one site are coupled ferromagnetically, (si(n) 
7r/4 < \a\ < 7r/2 the correlation function F{r) is always negative: 



F(r) 



Sj(n)) > 0. On the other hand, in the range 



(78) 



and all the correlation functions at one site are also negative. 

These ranges have two end points in common, a = ±7r/4, where Tc ~ 0. Whereas a = tt/A corresponds to the 
trivial partition of the system into independeat singlet pairs, the case a = — 7r/4 is more interesting. In this case 
one can calculate all spin correlations exactlytil. The correlations of spins located on neighboring 'sites' of lattice at 
a — — 7r/4 are antiferromagnetic, while all other correlators are zero. 

pihe Hamiltonian for model ( [73| ) has a very combersome form and for the cases a — —tt/A and a = 7r/2 was given 

iiE. 

Our results suggest that the spin correlation functions decay exponentially with a correlation length ~ 1 for an 
arbitrary parameter a. We also assume that the decay of the correlation function is of the exponential type for the 
14-parameter model as well, i.e., for any choice of site spinor ^■^t^'^P. This assumption is supported in special cases: 
1) the partition of the system into one-dimensional chains with exactly known exponentially decaying correlation 
functions; 2) the two-dimensional AKLT model, for which the exponential character of the decay of the correlation 
function has been rigorously proveciE3. Further evidence of the stated assumption lies in the numerical results obtained 
for various values of the parameter in the one-parameter model. 
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FIG. 13. Example ol a graph corresponding to the wave function (|8l|). 

D. Generalization of the model to other types of lattices 

The wave function (^5|), ( |60|) can be generalized to any type of lattice. The general principle of wave function 
construction for a system of spins 1/2 entails the following: 

1) Each bond on a given lattice has associated with it two indices running through the values 1 and 2, one at each 
end of the bond. 

2) Each bond has associated with it a metric spinor g\fj_ with the indices of the ends of this bond. 

3) Each site of the lattice (a site being interpreted here, of course, in the same sense as in Sec.IVB) with m outgoing 
bonds has associated with it an mth-rank spinor with the indices of the bonds adjacent to the site. 

4) The wave function is the product of all spinors at sites of the lattice and all metric spinors. 

It is obvious that each index in the formulated wave function is encountered twice, so that the wave function is 
scalar and, hence, singlet. 

The wave function so constructed describes a system in which each lattice site contains as many spins s — 1/2 as 
the number of bonds emanating from it. 

To completely define the wave function, it is necessary to determine the specific form of all site spinors. The 
coefficients that determine their form are then parameters of the model. 

The Hamiltonian of such a model is the sum of the cell Hamiltonians acting in the spin space of the subsystem 
formed by the spins at two mutually coupled sites: 

H = Y,m,. (79) 

Each cell Hamiltonian is the sum of the projectors with arbitrary positive coefficients onto all multiplets possible 
in the corresponding two-site subsystem except those present in the constructed wave function: 

H^., ^Y^^kPi:' ■ (80) 

k 

Then = and, accordingly, s) = 0. 

Consequently, ^ s is an exact ground-state wave function. 

We note that any two lattice sites can be joined by two, three, or more bonds, because this does not contradict the 
principle of construction of the wave function. Moreover, the general principle of construction of the wave function is 
valid not only for translationally symmetric lattices, but for any graph in general. As an example, let us consider the 
system shown in Fig. |l^. The wave function of this system has the form 

VP, = M/^i (1)VI/A.M1-1P1 {2)^P.^'.r^ (3)^M2r. (4)g;,,,,5^,^,5,,,,g,,,,5,,,, (81) 

and describes a system containing ten spins 1/2. 

If the given lattice has dangling bonds (as occurs for systems with open boundary conditions), the resulting wave 
function represents a spinor of rank equal to the number of loose ends. The ground state of this kind of system is 
therefore 2'-fold degenerate, where I is the number of loose ends. For an open one-dimensional chain, for example, 
the ground state corresponds to four functions — one singlet and three triplet components. For higher-dimensional 
lattices this degeneracy depends on the size of the lattice and increases exponentially as its boundaries grow. 
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V. ELECTRONIC MODELS 



In r e j Cent . years there has been increasing interest in studying models where at least the ground state can be found 
exactlyBEJ. The_most popular method for the construction of exact ground state is so-called optimal ground state 
approach (0GS)E3. In the OSG method the ground state of the system is simultaneously the ground state of each 
local interaction. In this Section we propose new ID and 2D models of interacting electrons with the exact ground 
state which are different from those constructed in the OSG method. The ground state wave function of our models 
is expressed in terms of the two-particle 'singlet bond' (SB) function located on sites i and j of the lattice: 

J'] = 4Ai - <i4T + ^ «T<i + 4A..i) io) ' (82) 

where cf^, Ci^p- are the Fermi operators and x is an arbitrary coefRcient. The SB function is the generalization of the 
RVB functionEj including ionic states. The presence of the ionic states is very important from the physical point of 
view because, as a rule, the bond functions contain definite amount of the ionic states as well. 

It is knawnj-a. iS^^-pf ID and 2D quantum spin models the exact ground state of which can be represented in the 
RVB fornrt3BE3^t3^Eil. It is natural to try to find electronic models with exact ground state at half-filling formed by 
SB functions in the same manner as for above mentioned spin models. The electronic models of these types include 
the correlated hopping of electrons as well as the spin interactions and pair hopping terms. 

The model with dimerization. 

As the first example we consider the ID electronic model with the two-fold degenerate ground state in the form of 
the simple product of SB dimers, similarly to the ground state of the well-known spin-i Majumdar- Ghosh modelQ. 
For the half-filling case the proposed ground state wave functions are: 

«'i = [l,2][3,4]...[iV-l,7V] (83) 

and 



*2 - [2,3][4,5]...[iV,l] (84) 

In order to find the Hamiltonian for which the wave functions ( |83| ) and ( |8^ ) are the exact ground state wave 
functions, we represent the Hamiltonian as a sum of local Hamiltonians hi defined on three neighboring sites (periodic 
boundary conditions are supposed): 

N 

H = Y,K (85) 
i=i 

The basis of three-site local Hamiltonians hi consists of 64 states, while only eight of them are present in 5'i and 
^2- These 8 states are 

[i,i + l](p,+2, [i + l,z + 2], (86) 

where ipi is one of four possible electronic states of i-th site: |0)j, It)^, U),;, |2)j. 

The local Hamiltonian hi for which all the functions (^6|) are the exact ground state wave functions can be written 
as the sum of the projectors onto other 56 states \xk) 

K=Y,^k\Xk){Xk\, (87) 

k 

where \k are arbitrary positive coefhcients. This means that the wave functions ^'i and ^'2 are the ground states 
of each local Hamiltonian with zero energy. Hence, VPi and 5*2 ye. the. 'nptimar ground state wave functions of the 
total Hamiltonian H with zero energy, similarly to the models irJlJ'oEZlllj. In general case, the local Hamiltonian hi 
is many-parametrical and depends on parameters Afe and x. We consider one of the simplest forms of hi including the 
correlated hopping of electrons of different types and spin interactions between nearest- and next-nearest neighbor 
sites: 

hi ~ 2 — X (ii,i+i + ii+i.i+2) 

+ (a;2-(l + x2)(l-n,+i)2)T,,,+2 
1 — a;2 

+ 8 {Si ■ Si-|_i + Si_|_i • S,;+2 + Si • Si+2), 
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FIG. 14. The lattice ol the Shastry - Sutherland model 



where 

^i^J = ^^4,o'^3,'y + Cj>c*,^)('^»,-<T - nj,-af (89) 

a 

and Si is the SU{2) spin operator. 

Each local Hamiltonian hi is a non- negatively defined operator at < 1. The following statements related to the 
Hamiltonian (|88| ) are valid: 

1. The functions ( |83| ) and ( p^ ) are the only two ground state wave functions of the Hamiltonian ( |88| ) at A^e = 
(iVe is the total number of electrons). They are not orthogonal, but their overlap is ~ at 3> 1. 

2. The ground state energy Eq^N^/N) is a symmetrical function with respect to the point Ni,/N = 1 and has a 
global minimum _Eo = at N^/N — 1. 

3. The translational symmetry of ( p8[ ) is spontaneously broken in the ground state leading to the dimcrization: 

{\ti,i+l — ti+is+2\) = 2 

The excited states of the model can not be calculated exactly but we expect that there has to be a gap, because 
the ground state is formed by the ultrashort-range SB functions. If it is the case, the function EQ{Ne/N) has a cusp 
at Ne/N = 1. 

Actually, this model is the fermion version of the Majumdar - Ghosh spin model. Moreover, it reduces to the 
Majumdar - Ghosh model at a; = and in the subspace with Ui = 1. 



For X — I the Hamiltonian (88) simplifies and takes the form: 

H = -2J2 E (90) 

3 3 

The 2D model. 

We can easily construct the 2D electronic model with the exact ground state which is analogous to the Shastry - 
Sutherland modelEj (Fig.jlJ). The Hamiltonian of this model is: 

H= ^^,3+Kk + hlu. (91) 

where the sum is over all triangles {i, j. A:}, one of which is shown on Fig.|lJ. So, each diagonal line belongs to the two 
different triangles. The local Hamiltonians h'^ f. acting on the diagonal of the triangle {ijj. A:}, and /ij^j, hi^k have the 
form (for the sake of simplicity we put a; = 1) 



Kk = -t,- fc - e^""^ T,- (92) 



It is easy to check that 
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All other states of the local Hamiltonian hij + hi^k + ^^fc have higher energies. Therefore, the ground state wave 
function in the half-filling case is the product of the SB functions located on the diagonals shown by dashed lines on 
Fig.jlJ. This model has non-degenerate singlet ground state with ultrashort-range correlations. 

The ladder model. 

Let us now consider electronic models with a more complicated ground state including different configuipajtions 
of short-rango-SB functions. The form of these ground states is similar to that for spin models proposed intS and 
generalized inEJ. In the ID case our model describes two- leg ladder model (Fig.^). Its ground state is a superposition 
of the SB functions where each pair of nearest neighbor rungs of the ladder is connected by one SB. 

The wave function of this ground state can be written in a form (pq): 



The functions tp^^ii) describes i-th rung of the ladder 



with 



/IT)fe\ 

\l)k 

|2>. , 

V |o), / 



It is easy to see that 
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(93) 



(94) 



(95) 



Therefore, the function 'i's is a singlet wave function depending on two parameters x and Ci/c2. Actually, this form 
of '^s is equivalent to the MP form with 4x4 matrices Ax,y{i) = gxf^ip^'^ (i). Moreover, at x = and ci/c2 = —1 the 
function 'i's reduces to the wave function of the well-known AKLT spin-1 model. 

In order to find the Hamiltonian for which the wave function (|9^) is the exact ground state wave function, it is 
necessary to consider what states are present on the two nearest rungs in the "^g. It turns out that there are only 16 
states from the total 256 ones in the product il^^^{i)gpvilj^''{i + 1). The local Hamiltonian hi acting on two nearest 
rungs i and i -f 1 can be written in the form of ( p7| ) with the projectors onto the 240 missing states. The total 
Hamiltonian is the sum of local ones ( p5| ) . The explicit form of this Hamiltonian is very cumbersome and, therefore, 
it is not given here. 

The correlation functions in the ground state ( p3| ) can be calculated exactly in the same manner as it was done for 
spin modelalj. It can be shown that all of correlations exponentially decay in the ground state. We expect also that 
this model has a gap. 

This method of cottstruction of the exact ground state can be generalized also to 2D and 3D lattices, as it was done 
in Sec. IV. FoUowingEil, one can rigorously prove that the ground state of these models is always a non-degenerate 
singlet. 

ID models with the giant spiral order. ■— ■ 

There is one more spin-i model with an exact ground state of the RVB typeL3. This is the model (^ describing the 
F-AF transition point. The exact singlet ground state can be expressed by the combinations of the RVB functions 
(i,j) distributed uniformly over the lattice points (||). The analog of the wave function (||) in the SB terms is: 



E 

i<j... 



(-1)^ [i,j][kj][m,n] 



(96) 



where P = k,l, . . .) is the permutation of numbers (1, 2, . . . N). It is interestiHjg to note that the singlet wave 
function ( p6| ) can be also written in the MP form but with an infinite size matricesEJ. 

In order to find the Hamiltonian for which the wave function (96) is the exact-ground state wave function, let us 
consider what states are present on the two nearest sites in the ^'q. It turns outO that there are only 9 states from 
the total 16 ones in (Bo). They are 



ITT), lii), ITi + iT), 120 - 02) 
ITT -TT) + 2; 120 + 02), IT 0-0 T), 
IT2-2T), ITO-Oi), IT 2-2 i) 



(97) 
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The local Hamiltonian /i^^i+i can be written as the sum of the projectors onto the 7 missing states (|87|). At \x\ > 1 
the most simple form of this total Hamiltonian is: 



N 



H — I 2^1, j+i — ^ 4 Si • Si+i 



4 .t2 

H ^■'yi • Vi+1 + 4 ■ 



(98) 



We use here rj operators: 



1- Ui 



which form another SU{2) algebraSEi, and iji • 772 is a scalar product of pseudotsspins rji and 772- We note that this 
Hamiltonian commutes with S^, but does not commute with r]^. It can be provedED that only three multiplets are the 
ground states of (p8|): the singlet state (p6|), the trivial ferromagnetic state S = N/2 and the state with 5* = NJ2 — 1. 

The norm and the correlators of the electronic model (|9^) in the singlet ground state are exactly calculatedliZl. For 
example, the norm of (pq) is: 



J AT 



N 



^ 1 + cosh(a;^) 
cos2(|) 



The correlators at ^ 1 are 



(nM^^ = o{ — 



Cr+l,a) = O ( — 



(S,Si+/> = icos 



(99) 



The correlator {r]fri^_^i) which determines the off-diagonal long-range order (ODLRO)@ vanishes in the thermody- 
namic limit. At the same time the spin-spin correlations have a spiral form, and the period of the spiral equals to the 
system size as in the spin model (|9|). 

Another electronic model can be obtained by making the canonical transformation c^^ — s- c^^ and c^j^ — s- Ci^i- As 
a result of this transformation, the SB function (82) becomes: 



and the wave function ( p6| ) changes to 

*o = XI {hj}{kA{m,n} ... 



(100) 



(101) 



The function (101) for |a;| > 1 is the exact ground state wave function of the transformed Hamiltonian: 

N 



1=1 

4 

^S,.S,+i+4- 



i ■ Vi+l 



-S^Sl^i + 1 



(102) 



where 



This Hamiltonian commutes with 77^ but does not commute with S^. Therefore, the eigenfunctions of the Hamiltonian 
( |102| ) can be describied by quantum numbers nand 77^ . For the cyclic model the state s wit h three different values of 
77 have zero energytZl [as it was for the model (Pq)]. They include one state with 77 = (101), all states with 77 — N/2: 
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N/2, ri' 



+ \N/2-7j' 



|0) 



(103) 

0), the ground state of the 



and th e sta tes with t] = N/2 — 1. Therefore, for the case of one electron per site (ry^ 
model (102) is three- fold degenerate. 

The correlation functions in the ground states with 77 — N/2 and 77 — N/2 — 1 for the half- filling case coincide with 
each other and at iV ^ 1 they are: 



(SiSi+i) 



1 



(104) 



The e xiste nce of the ODLRO immediately follows from the latter equations. The correlation functions in the ground 
state (101) have similar forms as in Eqs.(|99|): 



{eta = O 



N 



= o 



1 

iV2 



■ cos 



2t:1 



(105) 



The giant spiraL— ordering in the las t eq uation implies the existence of the ODLRO and, therefore, the 
superconductivityo in the ground state 



Similarly to the original spin model 



the last two electronic models (|98|),(102) describe the transition points 
on the phase diagram between the phases with and without a long-range order [ferromagnetic for the model ( |9^ ) 
and off-diagonal for the model (102)]. Therefore, we suggest the formation of the ground state with giant spiral 
order (ferromagnetic or off-diagonal) as a probable scenario of the subsequent destruction of the ferromagnetism and 
superconductivity. 



VI. CONCLUSION 



We have considered the class of the ID and 2D spin and electronic models with the exact ground states. 

The one of these models is the spin-^ ladder with competing interactions of the ferro- and antiferromagnetic types 
at the F-AF transition line. The exact singlet ground state wave function on this line is found in the special form 
expressed in terms of auxiliary Bose-operators. The spin correlators in the singlet state show double-spiral ordering 
with the period of spirals equal to the system size. 

In general case the proposed form of the wave function corresponds to the MP form but with matrices of infinite size. 
However, for special values of parameters of the model it can be reduced to the standard MP form. In particular, we 
consider spin-i ladder with nondegenerate antiferromagnetic ground state for which the ground state wave function 
is the MP one with 2x2 matrices. This model has some properties of ID AKLT model and reduces to it in definite 
limiting case. 

The ground state wave function of the spin ladders can be represented in an alternative form as a product of 
second-rank spinors associated with the lattice sites and the metric spinors corresponding to bonds between nearest 
neighbor sites. Two-dimensional spin-i model is constructed with exact ground state wave function of this type. The 
ground state of this model is a nondegenerate singlet with exponentially decay of spin correlators. We believe the 
model has a gap in the spectrum of excitations. 

We propose new models of interacting electrons with the exact ground state formed by the singlet bond functions 
in the same manner as for some spin models. In particular, we considered the models describing the boundary points 
on the phase diagram between the phases with and without long-range order (ferromagnetic or off-diagonal). 

In conclusion we note that the construction of considered models is based on following property. Their Hamiltonians 
are the sums of the Hamiltonians that are local and non-commuting with each other. In the same time the ground 
state wave function of the total Hamiltonian is the ground state one for each for them. It is clear that these models 
are rather special. Nevertheless, the study of them is useful for understanding properties of the real frustrated spin 
systems and strongly correlated electronic models. 

This work was supported by the Russian Foundation for Basic Research (grants no. 00-03-32981, 99-03-3280 and 
no.00-15-97334). 
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